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Abstract 

In this paper we calculate explicitly the classical secular precessions of the node Q and the perigee u> of 
an Earth artificial satellite induced by the even zonal harmonics of the static part of the geopotential up 
to degree I = 20. Subsequently, their systematic errors induced by the mismodelling in the even zonal 
spherical harmonics coefficients J\ are compared to the general relativistic secular gravitomagnetic and 
gravitoelectric precessions of the node and the perigee of the existing laser-ranged geodetic satellites 
and of the proposed LARES. The impact of the future terrestrial gravity models from CHAMP and 
GRACE missions is discussed as well. Preliminary estimates with the recently released EIGEN-1S 
gravity model including the first CHAMP data are presented. 
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1 Introduction 

Recently, great efforts have been devoted to the investigation of the possibility of measuring 
some tiny general relativistic effects in the gravitational field of the Earth by analyzing the 
laser-ranged data to some existing or proposed geodetic laser-tracked (SLR) satellites. 

The most famous experiment is that peformed with LAGEOS and LAGEOS II [Ciufolini 
et al, 1998] and aimed to the detection of the gravitomagnetic Lense-Thirring drag of inertial 
frames [Lense and Thirring, 1918; Ciufolini and Wheeler, 1995] in the gravitational field of 
the Earth. The analysis of the orbits of the LAGEOS satellites could allow also for an alter- 
native measurement of the gravitoelectric perigee advance [Ciufolini and Wheeler, 1995] to be 
performed in the gravitational field of the Earth [Iorio, 2002; Iorio et al., 2002a]. Moreover, 
the possibility of including also the data from other existing SLR satellites in these analysis is 
currently investigated [Iorio, 2002]. The proposed LAGEOS-LARES mission [Ciufolini, 1986], 
whose original configuration is currently being reanalyzed [Iorio et al, 2002b] in view of the 
inclusion of more orbital elements of various SLR satellites in the observable to be adopted, 
should be of great significance for both gravitomagnetic and gravitoelectric tests [Iorio et al, 
2002a; Iorio et al, 2002b]. Satellite laser ranging could be the natural candidate also for the 
implementation of a space-based experiment aimed to the detection of the so called gravit- 
omagnetic clock effect [Mashhoon et al, 1999; Iorio et al, 2002c], which is sensitive to the 
direction of motion of two counter-orbiting satellites along identical orbits in the gravitational 
field of a central rotating mass. 

In all such performed or proposed experiments it is of the utmost importance to reliably 
assess the error budget. Indeed, the terrestrial space environment is rich of competing classical 
perturbing forces of gravitational and non-gravitational origin which in many cases are far 
larger than the general relativistic effects to be investigated. In particular, it is the impact of 
the systematic errors induced by the mismodelling in such various classical perturbations which 
is relevant in determining the total realistic accuracy of an experiment like those previously 
mentioned. 

The general relativistic effects of interest here are linear trends affecting the perigee to and 
the node Q of the orbit of a satellite and amounting to 10 1 -10 3 milliarcseconds per year (mas/y 
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in the following) for the gravitomagnetic and the gravitoelectric effects, respectively. 

In this context the most important source of systematic error is represented by the sec- 
ular classical precessions of the node and the perigee induced by the mismodelled even (I = 
2n, n = 1,2, 3, ...) zonal (m = 0) harmonics 5J2, SJ4, 5Jq, ... of the multipolar expansion of the 
terrestrial gravitational field, called geopotential. Indeed, while the time-varying orbital tidal 
perturbations [Iorio, 2001; Iorio and Pavlis, 2001; Pavlis and Iorio, 2002] and non-gravitational 
orbital perturbations [Lucchesi, 2001; 2002], according to their periods P and to the adopted 
observational time span T Q b s , can be viewed as empirically fitted quantity and can be removed 
from the signal, this is not the case of the classical even zonal secular precessions and of certain 
subtle non-gravitational secular effects of thermal origin [Lucchesi, 2002]. Their mismodelled 
linear trends act as superimposed effects which may alias the recovery of the genuine general 
relativistic features. Such disturbing trends cannot be removed from the signal without can- 
celling also the general relativistic signature, so that one can only assess as more accurately as 
possible their impact on the measurement. Then, the systematic error induced by the mismod- 
elled part of the geopotential can be viewed as a sort of unavoidable, lower bound of the total 
systematic error. 

In this paper we calculate explicitly, up to I = 20, the expressions of the coefficients of 
the classical secular precessions on the node and the perigee due to the geopotential (section 
2). Their explicit, analytic form, although rather cumbersome, may turn out to be useful in 
designing suitably alternative relativistic observables which are not sensitive, at least in part, to 
such classical aliasing effects [Iorio and Lichtenegger, 2002; Iorio and Lucchesi, 2002]. In section 
3 we work out the numerical values of the mismodelled precessions for the existing SLR geodetic 
satellites and of the proposed LARES and compare them to the general relativistic effects. The 
errors for the spherical harmonics coefficients are those of EGM96 gravity model [Lemoine et 
ai, 1998] and of the recent EIGEN-1S (see |http: / / op.gfz-potsdam.de/ champ/results/p which 
includes the first data from CHAMP mission^]. These estimates can be useful in assessing 
the systematic errors in various possible observables, built up with such Keplerian orbital 
elements, which are sensitive to some relativistic effects. E.g., someone could look at the 

1 It should be pointed out that the values employed in the following for the errors 8J1 by EIGEN-1S are the 
formal, uncalibrated standard deviations. Moreover, EIGEN-1S is based only on the satellites tracking data. 
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Table 1: Orbital parameters of the existing spherical passive geodetic laser-ranged satel- 
lites and of LARES. Aj=Ajisai, Stl=Stella, Str=Starlette, WS=WESTPAC1, E1=ETAL0N1, 
E2=ETALON2, L1=LAGE0S, L2=LAGE0S II, LR=LARES. a is in km, i in deg and n in 
s" 1 . 





Aj 


Stl 


Str 


WS 


El 


E2 


LI 


L2 


LR 


a 


7,870 


7,193 


7,331 


7,213 


25,498 


25,498 


12,270 


12,163 


12,270 


e 


0.001 





0.0204 





0.00061 


0.00066 


0.0045 


0.014 


0.04 


i 


50 


98.6 


49.8 


98 


64.9 


65.5 


110 


52.65 


70 


n 


0.0009 


0.001 


0.001 


0.001 


0.00015 


0.00015 


0.00046 


0.00047 


0.00046 



relativistic perturbations of the radial, along-track and cross-track components of the position 
and velocity vectors of a satellite; they can be obtained from suitable combinations of the 
Keplerian orbital elements. The same holds also for the range and range-rate perturbations in 
intersatellite tracking missions like GRACE [Cheng, 2002]. In section 4 we review the strategy of 
combining the orbital residuals of the nodes and the perigees of different laser-ranged satellites 
which allows for a reduction of the impact of the geopotential's error and yield preliminary 
estimates of the errors affecting two gravitomagnetic combinations based on the first results 
from EIGEN-1S. They are useful in order to get an insight of the improvements which will 
take place when the full new gravity models will become available. Section 5 is devoted to the 
conclusions. 

In Tab. 1 we quote the orbital parameters of the existing spherical passive geodetic laser- 
ranged satellites Ajisai, Stella, Starlette, WESTPAC1, ETALON1, ETALON2, LAGEOS, LA- 
GEOS II and of the proposed LARES. In it a is the semimajor axis, e is the eccentricity, i is the 
inclination and n = V GMa~ 3 , where G is the Newtonian gravitational constant and M is the 
mass of the central body, is the Keplerian mean motion. It is worth noting that the perigees 
of many of them, except for Starlette, cannot be employed for any relativistic tests due to the 
notable smallness of their eccentricities. 
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2 The orbital classical precessions 

Here we show the explicitly calculated coefficients 

nrSeven zonal 

O _ u "class /i\ 

^ = d(J 2n ) (1) 

and 

fi. -even zonal 

^ = d(J 2n ) (2) 

of the satellites' classical nodal and apsidal precessions due to the even zonal harmonics of the 
geopotential up to / = 20. The classical precessions of the node and the perigee due to the even 
zonal harmonics of geopotential can be written as 



^s™ 1 = E^«X J 2n, ( 3 ) 

71=1 



CO. 



even zonal 
class 



= U.2n X J2n- (4) 



n=l 



As we shall see later, the coefficients Q.2n and u,2n depend only on the orbital parameters of the 
satellites. Recall that J\ = —Ciq, I = 2n, n = 1,2,3... where the unnormalized adimensional 
Stokes coefficients C\ m of degree / and order m can be obtained from the normalized C/ m with 

Clm = Ni m Cl m . (5) 

In it 

\ (2l + l)(2-6 0m )(l-m)\ ] * 

Nlm= [ (T+m)\ J • (6) 

The general expressions of the classical rates of the near Earth satellites' Keplerian orbital 
elements due to the geopotential d c i ass , e c i ass , M ass , ^ciass, ^ciassi -^ c i ass , and of the inclination 
functions Fi mp (i) and of the eccentricity functions Gi pq (e) can be found in [Kaula, 1966]. The 
coefficients Q,2n and u_2n are of crucial importance in the evaluation of the systematic error due 
to the mismodelled even zonal harmonics of the geopotential; moreover, they enter the combined 
residuals' coefficients q about which we speak in section 4. Since the general relativistic effects 
investigated are secular perturbations, we have considered only the perturbations averaged over 
one satellite' s orbital period. This has been accomplished with the condition / — 2p + q = 
which allows for canceling out the rate of the mean anomaly M.. Since the eccentricity functions 
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Gi pq are proportional to e' q \ for a given value of / we have considered only those values of p 
which fulfil the condition / — 2p + q = with q = 0, i.e. p — |. This implies that in the 
summations 

E E (7) 



i +oo J/O 

E^«* E ^? (8) 



and 

l +oo Q 

E F w P E ~ dp 

p=0 q=-oo ue 

involved in the expressions of the classical rates we have considered only F m i_ and G^ . More- 
over, in working out the G^q we have neglected the terms of order 0(e k ) with k > 2. 

2.1 The nodal coefficients 

The nodal coefficients, proportional to 

' E g 1m j:^, (9) 
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2.2 The perigee coefficients 

The coefficients of the classical perigee precession are much more involved because they are 
proportional to 

A E GtoZ^ + V—il E ^*E*W (20) 



We can pose dj. 2n = u a 2n + u b 2n . 

The first set is given by (R is the Earth's mean equatorial radius) 



u a 2 



/'i a 



3 

2 \ a 



R\ 2 cos 2 1 



UJ A = OJ 



.2 



a;_ 2 



1 -e 2 ) 2 ' 
5 /.RXMI + Ie 2 ) 



a / 



=2X2 



(7 sin 2 i- 4) 



35 /R\ 4 (l + 5e 2 



a; s = a> 



.2 



(1 - e 2 ) 4 

21 
2 



33 a 9 o 

— sin ? sin i 

8 2 



105 (R\ & (1 + fe 2 ) /715 . 



16 V a7 (1-e 2 ) 6 V64 



sin 6 i- 



(21 

(22) 

(23) 



w .io = 



143 . , . 33 




■sin i H 


8 


4 




"1,155 


,'\ a 






128 W 



'1-e 



2X8 



4,199 . o . 1,105 , 

sin % sin i 

128 16 



(24) 



+ 



195 



L0 12 = U 



.2 



sin i — 13 sin i + 1 



3,003 fR\ 10 (1 + fe 2 ) /52,003 . 



_ 256 W (1-e 2 ) 10 V 512 



in. 33^915 . q. 
sin i sin % 



128 



(25) 



8,075 , 425 4 75 2 
+ ^ sin 2 — sin % + — sm z — 1 



W .14 — W .2 



15, 015 fR\ 12 (1 + f e 2 ) /334, 305 _ ; 



1,024 \aj (1-e 2 ) 12 V 1,024 



12 260,015 in 

sm i — - — sm u i+ 

256 



(26) 



156,009 o 11,305 e 1,615 4 
+ — 7777 — sm i — — sm % H — sm %— 



128 



16 



8 



16 



51 



sin 2 z 



36,465 (R 



2,048 \a 



ii 



1 + 



l-e 



2U4 



17,678,835 . 14 . 3,991,995 . 12 . 

sin z — , — sm z 



16,384 



1,024 



(27) 



2,890,755 . in . 535,325 . , 107,065 . fi 
+ ' sm 10 z --— — sm 8 z -\ ^ sm 6 z 



512 



128 



64 



2,793 4 133 2 

sm z H — sm i — 1 

8 4 



w .i8 = 



692,835 (R\ 16 (1 + 68e 2 ) /119, 409, 675 . 



32,768 Va/ (l-e 2 ) 16 \ 32,768 



sm 16 z 



(28) 



30,705,345 . u . 6,513,255 . 12 . 1,470,735 . 10 . 

sm 14 1 + ' sm 12 i - — -\ sm 10 z + 



2,048 



256 



64 



760,725 ,. 28,175 . fi . 1,127 . 4 . 

+ sin z sin ? H sm z- 

64 8 2 



42 sin 2 z + 1 



20 



.2 



1,616,615 /R\ 18 (1 + i|ie 2 ) /l, 641, 030, 105 18 



65,536 \aj (l-e 2 ) 18 \ 131,072 



(29) 



1,893,496,275 . lfi . 460,580,175 . u . 30,705,345 . 12 . 

■an H . sm z sm z 



32, 768 



4,096 



256 



19,539,765 . in . 1,890,945 . 8 . 108,675 . , 
+ — hr^ sin 10 z - -2— -' sin 8 z + — ^ — sin 6 z- 



256 



64 



16 



1,725 . 4 . 207 . 2 . 1 

— sm z H sm z — 1 

2 4 

The second set is given by (R is the Earth's mean equatorial radius) 
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3 The mismodelled classical precessions 



The results obtained in the previous section can be used in working out explicitly the contribu- 
tions of the mismodelled classical nodal and apsidal precessions up to degree I = 20 of the exist- 
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ing spherical passive laser-ranged geodetic satellites and of the proposed LARES. They are of the 
form 5ti( 2n ) = ^.2n x dJ2n, n = 1,2, ...10 and 5to^ n ) = cl>. 2 n x <5J 2 „, n— 1,2, ...10. The coefficients 
Cl,2n and cj.2n are worked out in section 2 and the values employed for 5J 2n = — y/An + lx<5C 2n o, 
n = 1, 2, ...10 are those quoted in the adopted Earth's gravity model. 

Table 2: Mismodelled classical nodal precessions SQ^n) and predicted Lense-Thirring nodal 
precessions ^lt of the existing geodetic laser-ranged satellites and of LARES. Ll=LAGEOS, 
L2=LAGEOS II, LR=LARES, Aj=Ajisai, Stl=SteUa, Str=Starlette, WS=WESTPAC1, 
El=ETALONl, E2=ETALON2. All the values are in mas/y. For the ETALON satellites, 
when the values are less than 10" 4 mas/y a - has been inserted. EGM96 gravity model has 
been adopted. 
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From Tab. 2 it is interesting to note that for the satellites orbiting at lower altitudes than 
the LAGEOS satellites the impact of the mismodelled part of the geopotential does not reduce 
to the first two or three even zonal harmonics. This feature is very important in calculating the 
error budget, especially if the nodes of low-orbiting satellites are to be considered. Moreover, 
while for the LAGEOS family a calculation up to I = 20 is rather adequate, this is not the case 
for the other satellites for which the even zonal harmonics of degree I > 20 should be considered 
as well. In regard to this topic, the choice of the Earth gravity model becomes crucial because 
EGM96, for example, does not seem to be particularly reliable at degrees higher than 20. The 
same considerations hold also for the perigee whose mismodelled classical precessions are quoted 
in Tab. 3. We have considered only LAGEOS II, Starlette and the LARES due to the extreme 
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smallness of the eccentricity of the other satellites. From both Tab. 2 and Tab. 3 the relevant 

Table 3: Mismodelled classical perigee precessions 5Cj{2n) and predicted Lense-Thirring and 
gravitoelectric perigee precessions cOlt an d c^ge of the existing spherical passive geodetic laser- 
ranged satellites and of LARES. Ll=LAGEOS, L2=LAGEOS II, LR=LARES, Aj=Ajisai, 
Stl=SteUa, Str=Starlette, WS=WESTPAC1, El=ETALONl, E2=ETALON2. All the values 
are in mas/y. EGM96 gravity model has been adopted 



2n 


LI L2 


LR 


Aj 


Stl Str WS El E2 


2 


-42.3 


20.3 




-320.7 


4 


-122.7 


-17.6 




-1,924.4 


6 


-18.2 


-49.2 




429.1 


8 


43.1 


-42.6 




6,355.8 


10 


19.5 


-18 




2,805.1 


12 


-5.3 


-3 




-10,862.2 


14 


-6.2 


2 




-10,774.7 


16 


-0.2 


1.3 




8,395.8 


18 


0.4 


0.4 




9,086.4 


20 


0.1 


0.08 




-3,043.3 




-57.5 


-31.6 




68.5 




3,348 


3,278.6 




11,804.7 



impact of the first two or three even zonal harmonics, at least for the LAGEOS satellites, is 
quite apparent. The situation for the lower orbiting satellites is far more unfavorable. 

In Tab. 4 and Tab. 5 we repeat the analysis with the very recently released EIGEN-1S 
gravity model which includes the first data from CHAMP. 

From an inspection of Tab. 2-Tab. 5 it turns out very clearly that, if the orbital elements 
of satellites other than those of LAGEOS family are to be considered, the observable which 
would account for them must cope with the problem of reducing the impact also of the degrees 
higher than 4. 

4 The systematic zonal error 

A possible strategy for reducing the impact of the geopotential's error consists of suitable 
combinations of the orbital residuals of the rates of the nodes and the perigees of different SLR 
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Table 4: Mismodelled classical nodal precessions Sfl^n) and predicted Lense-Thirring nodal 
precessions f^T of the existing geodetic laser-ranged satellites and of LARES. Ll=LAGEOS, 
L2=LAGEOS II, LR=LARES, Aj=Ajisai, Stl=Stella, Str=Starlette, WS=WESTPAC1, 
El=ETALONl, E2=ETALON2. The errors 5J 2n are those of the preliminary EIGEN-1S Earth 
gravity model from 88 days of CHAMP data. All the values are in mas/y. For the ETALON 
satellites, since all the values are less than 10 _1 mas/y a - has been inserted. 



2n 


LI 


L2 


LR 


Aj 


Stl 


Str 


WS 


El 


E2 


2 


-3.9 


7.2 


3.9 


35 


-11.1 


45.1 


-10.2 






4 


-7.2 


2.6 


7.3 


7.7 


-77.7 


8.9 


-71.7 






6 


-3.8 


-5.8 


3.8 


-181.5 


-204.5 


-313.4 


-190.1 






8 


-0.4 


-2.3 


0.4 


-83.8 


-340.3 


-154.3 


-320.3 






10 


0.4 


0.6 


-0.4 


168.5 


-379.6 


395.6 


-364.6 






12 


0.3 


0.4 


-0.3 


125.9 


-417 


318.2 


-413.5 






14 


0.1 


-0.001 


-0.1 


-122.1 


-441.9 


-400.8 


-461.3 






16 


0.02 


-0.07 


-0.03 


-152 


-354.1 


-525.6 


-404.7 






18 


-0.003 


-0.01 


0.003 


49.2 


-214 


242.2 


-299.4 






20 


-0.005 


0.007 


0.006 


131.3 


-42.2 


624.3 


-163.9 








30.7 


31.6 


30.8 


116.7 


152.8 


144.4 


151.5 


3.4 


3.4 



satellites [Ciufolini, 1996; Iorio, 2002]. Such combinations can be written in the form 

N 

(kfi = X GRf i GR , (42) 

8=1 

in which the coefficients Cj are, in general, suitably built up with the orbital parameters of 
the satellites entering the combinations, the /, are the residuals of the rates of the nodes and 
the perigees of the satellites entering the combination, Agr is the slope, in mas/y, of the 
general relativistic trend of interest and /xqr is the solve-for parameter, to be determined by 
means of usual least-square procedures, which accounts for the general relativistic effect. For 
example, in the case of the Lense-Thirring-LAGEOS experiment [Ciufolini, 1996] Alt = 60.2 
mas/y, while for the gravitoelectric perigee advance [Iorio, 2002] Age = 3,348 mas/y. More 
precisely, the combinations of eq. ( |4"2]) are obtained in the following way. The equations for the 
residuals of the rates of the N chosen orbital elements are written down, so to obtain a non 
homogeneous algebraic linear system of iV equations in A unknowns. They are ^gr and the 
first A — 1 mismodelled spherical harmonics coefficients 8 Ji in terms of which the residual rates 
are expressed. The coefficients q and, consequently, Agr are obtained by solving for /xgr the 
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Table 5: Mismodelled classical perigee precessions Scu^n) and predicted Lense-Thirring and 
gravitoelectric perigee precessions o>lt and uqe of the existing spherical passive geodetic laser- 
ranged satellites and of LARES. L1=LAGE0S, L2=LAGE0S II, LR=LARES, Aj=Ajisai, 
Stl=Stella, Str=Starlette, WS=WESTPAC1, E1=ETAL0N1, E2=ETALON2. All the val- 
ues are in mas/y. The errors SJ2 n are those of the preliminary EIGEN-1S Earth gravity model 
from 88 days of CHAMP data. 



2n 


LI L2 


LR Aj 


Stl Str WS El E2 
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-4.9 


2.4 


-37.8 


4 


-18.3 


-2.6 


-288.2 


6 


-4 


-11 


96.2 


8 


9.8 


-9.7 


1,454.6 


10 


4 


-3.6 


574 


12 


-1 


-0.6 


-2,112.2 


14 


-1.4 


0.4 


-2,425.3 


16 


-0.07 


0.3 


2,468.3 


18 


0.2 


0.1 


3,794.6 


20 


0.05 


0.04 


-1,614.6 




-57.5 


-31.6 


68.5 




3,348 


3,278.6 


11,804.7 



system of equations. So, the coefficients q are calculated in order to cancel out the contributions 
of the first N — 1 even zonal mismodelled harmonics which, as we have seen in the previous 
section, represent the major source of uncertainty in the Lense-Thirring and gravitoelectric 
precessions [Ciufolini, 1996; Iorio, 2002]. The coefficients c, can be either constant f\or depend 
on the orbital elements of the satellites entering the combinations through the coefficients Q,2n 
and u_2n worked out in section 2. 

Now we expose how to calculate the systematic error due to the mismodelled even zonal 
harmonics of the geopotential for the combinations involving the residuals of the nodes and the 
perigees of various satellites. 

In general, if we have an observable q which is a function q = q(xj), j = 1,2. ..M of M 

2 In general, the coefficient of the first orbital element entering a given combination is equal to 1, as for the 
combinations in [Ciufolini, 1996; Iorio, 2002]. 
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correlated parameters Xj the error in it is given by 

Sq 



M 



1 > 2 e (^Mf^W 2 



dx h 



dxi 



iik 



(43) 



in which crj = Cjj and a 2 ^ = where {Chk} is the square matrix of covariance of the 
parameters Xj. 

In our case the observable q is any residuals' combination 



N 



1 = E c ^( x i)' j = 1,2. ..10, 



(44) 



i=i 



where Xj, j = 1,2. ..10 are the even zonal geopotential's coefficients J2, J^...J%q. Since 



dq 



^ dfi 

E C »TT' J = l,2...10, 



(45) 



by putting eq. fl4"5| ) in eq. (^) one obtains, in mas/y 



Sq 



10 / N 



3=1 \i=l h^k=l \i=l 



' JV 

E- 



(46) 



h^fc=l \i=l ^""ft/ \j =1 ^fc, 

The percent error, for a given general relativistic trend and for a given combination, is obtained 
by taking the ratio of eq. (|46| ) to the slope in mas/y of the general relativistic trend for the 
residual combination considered. 

The validity of eq. ( flq) has been checked by calculating with it and the covariance matrix of 
EGM96 gravity model the systematic error due to the even zonal harmonics of the geopotential 
of the gravitomagnetic LAGEOS experiment; indeed the result 



5/lLT = 12.9% // LT 



(47) 



claimed in [Ciufolini et al., 1998] has been obtained again. For the EGM96-induced systematic 
error due to the even zonal harmonics of the geopotential of alternative proposed gravitomag- 
netic and gravitoelectric experiments, see [Iorio, 2002; Iorio et al., 2002a]. 

A very important point to stress is that the forthcoming new data on the Earth's gravita- 
tional field by CHAMP, which has been launched in July 2000, and GRACE, which has been 
launched in March 2002, will have a great impact on the reduction of the systematic error due 
to the mismodelled part of geopotential. 
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In order to get a preliminary insight of what the improvement due to the new Earth gravity 
models might be, let us consider, for example, the usual observable by Ciufolini [Ciufolini, 1996] 
for the detection of the Lense-Thirring drag 

5Q L1 + 0.295 x 5Q L2 - 0.35 x 5u L2 ~ 60.2/x LT . (48) 

The root-sum-square error due to geopotential, according to the diagonal part only of the 
covariance matrix of EGM96 model, amounts to 46.5% 0; according to the diagonal part only 
of the covariance matrix of EIGEN-1S (at present, its full covariance matrix is not yet publicly 
available), it reduces to 10.5%. 

If we consider the combination proposed in [Iorio, 2002] 

5(l hl + 0.444 x 5tt L2 - 0.027 x 5(l Ai - 0.341 x 5u h2 ~ 61.2/x LT , (49) 

in this case the root-sum-square error due to geopotential, according to the diagonal part only 
of the covariance matrix of EGM96 model, amounts to 64.2% Q according to the diagonal part 
only of the covariance matrix of EIGEN-1S, it reduces to 15.5%. 

It should be noted that, according to [Ries et al, 1998], it would not be entirely correct to 
automatically extend the validity of the full covariance matrix of EGM96, which is based on a 
multi-year average that spans the 1970, 1980 and early 1990 decades, to any particular time 
span like that, e.g., of the LAGEOS-LAGEOS II Lense-Thirring analysis which extends from 
the middle to the end of the 1990 decade. Indeed, there would not be assurance that the errors 
in the even zonal harmonics of the geopotential during the time of the LAGEOS-LAGEOS II 
Lense-Thirring experiment remained correlated exactly as in the EGM96 covariance matrix, 
in view of the various secular, seasonal and stochastic variations that we know occur in the 
terrestrial gravitational field and that have been neglected in the EGM96 solution. Then, the 
use of the diagonal part only of the covariance matrix of EGM96 should yield more conservative 
results. However, since it turns out that such seasonal effects would mainly affect just the first 
even zonal harmonic coefficients of the geopotential, the uncertainty related to them should be 

3 It reduces to 12.9% by considering also the correlation among the spherical harmonics coefficients, according 
to EGM96, as stated before. 

4 It reduces to 10.8% by considering also the correlation among the spherical harmonics coefficients, according 
to EGM96 [Iorio, 2002]. 
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very small for residual combinations which, by construction, cancel out just the first even zonal 
harmonic coefficients of the geopotential. On the other hand, if we cancel out as many even 
zonal harmonics as possible, the uncertainties in the evaluation of the systematic error based 
on the remaining correlated even zonal harmonics of higher degree should be greatly reduced, 
irrespectively of the chosen time span. This would have a relevant importance, e.g., for those 
even zonal harmonics like J 6 and J 8 whose favorable correlation in the covariance matrix of 
EGM96 seems to be the source of the perhaps optimistic evaluation of the systematic error 
due to the even zonal harmonics of the geopotential in the case of the LAGEOS-LAGEOS II 
Lense-Thirring experiment [Ries et ai, 1998]. The LAGEOS-LAGEOS II-LARES proposed 
combination of [Iorio et ai, 2002b] would cancel out, apart from 8J 2 and SJ4, just SJq and SJg. 

5 Conclusions 

The systematic error induced by the mismodelled static part of the geopotential is the ma- 
jor source of uncertainty in many proposed or performed tests of General Relativity in the 
gravitational field of the Earth via Satellite Laser Ranging. 

In this paper we have explicitly calculated the expressions of the coefficients of the classical 
secular precessions of the node and the perigee induced by the even zonal harmonics of the 
geopotential up to degree I = 20. The explicit expressions of the classical precessions may be 
useful, e.g., in getting insights for designing suitably new relativistic observables which are not 
too sensitive to such disturbing effects. 

Subsequently, we have compared the mismodelled precessions, according to EGM96 gravity 
model and EIGEN-1S preliminary gravity model which includes 88 days of data from CHAMP, 
to the general relativistic gravitomagnetic and gravitoelectric secular trends affecting the same 
orbital elements of the existing or proposed laser-ranged geodetic satellites. Since such satellites 
are the natural candidates for a number of relativistic tests in the gravitational field of the Earth, 
the presented calculations are useful in order to get an idea of the level of aliasing induced by 
the geopotential on the relativistic signatures of possible observables which may be built up 
with the orbital elements of such satellites. Of course, the obtained results could turn out 
to be useful also for other nonrelativistic investigations. While for the LAGEOS satellites a 
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calculation up to I — 20 is well adequate, for the other satellites orbiting at lower altitudes also 
the other harmonics of higher degrees should be carefully considered. The need for reducing 
the impact of the mismodelled classical precessions on the relativistic signals is quite apparent. 

Finally, we have shown how to calculate explicitly the static gravitational systematic error on 
suitably designed combinations involving the orbital residuals of different satellites. Preliminary 
estimates of the errors affecting such combinations with EGM96 and EIGEN-1S suggest that 
the future, more accurate terrestrial global gravity models from CHAMP and GRACE missions 
will have a notable impact on the improvement of, among other things, the precision of many 
general relativistic tests. 
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